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I begin to look at the growth of weak density inhomogeneities of nonrelativistic matter, in bimetric- 
MOND (BIMOND) cosmology. Far from making an exhaustive study, I concentrate on one attractive 
cosmological scenario, which employs matter-twin-matter-symmetric versions of BIMOND, and, 
furthermore, assumes that, on average, the universe is symmetrically populated in the two sectors. 
MOND effects are totally absent in an exactly symmetric universe, apart from the significant possible 
appearance of a cosmological constant, A ~ (ao/c)^. MOND effects-local and cosmological-do enter 
when density inhomogeneities that differ in the two sectors appear and develop. MOND later 
takes its standard form in systems that are islands dominated by pure matter, as are presumably 
the well scrutinized systems such as galaxies. I derive the nonrelativistic (weak-field-slow-motion) 
equations governing small-scale fluctuation growth. The equations split into two uncoupled systems, 
one for the sum, the other for the difference, of the fluctuations in the two sectors. The former is 
governed strictly by Newtonian dynamics, and describes standard growth of fluctuations. The 
latter is governed by MOND dynamics, which entails stronger gravity, and nonlinearity even for 
the smallest of perturbations. These cause the difference to grow faster than the sum, leading to 
anticorrelated perturbations, conducing to matter-twin-matter segregation (which continues for high 
overdensities) . The nonlinearity also causes interaction between nested perturbations on different 
scales. Because matter and twin matter (TM) repel each other in the MOND regime, matter 
inhomogeneities grow not only by their own self gravity, but also through shepherding by flanking 
TM overdensities (and vice versa). The relative importance of gravity and pressure in the MOND 
system (analog of the Jeans criterion), depends also on the strength of the perturbation. MOND 
gravity, which scales as the square root of the difference perturbation, holds sway over pressure for 
any mass, for weak enough perturbations. The development of structure in the universe, in either 
sector, thus depends crucially on two initial fluctuation spectra: that of matter alone and that of 
the matter-TM difference. I also discuss the back reaction on cosmology of BIMOND effects that 
appear as "phantom matter" (interpreted by some as "dark matter" ) , resulting from inhomogeneity 
differences between the two sectors. 

PACS numbers: 95.35.-|-d 



I. INTRODUCTION 

A full fledged treatment of cosmology and of structure formation in the MOND paradigm is still wanting. Structure 
formation on an expanding cosmological background, in nonrelativistic (NR) versions of MOND, has been studied 
with A^-body calculations, based on various heuristic extensions of MOND theories designed to work for isolated NR 
systems, e.g., in [1-5]. Clearly, a definite treatment of the problem has to be based on a fully relativistic theory, 
encompassing both cosmology and structure formation. Much work along such lines has been done in the framework 
of TeVeS-a relativistic MOND theory propounded by Bekenstein [6] (Sanders's stratified framework [7] is a precursor, 
and [8] is a recent review of TeVeS and its cosmology). For example, [9-14] considered various aspects of cosmology 
at large, the cosmic microwave background (CMB), and linear stages of structure formation in TeVeS-like theories. 
Although very encouraging, these results do not yet provide a fully satisfactory description of the phenomena at hand. 

A new relativistic theory for MOND (called BIMOND, for "bimetric MOND") was proposed recently based on a 
bimetric structure of space-time [15]. In this theory, gravity is described by two metrics instead of the one underlying 
general relativity (GR). Ordinary matter, to whose sector belongs our direct ken, couples, as in GR, only to one 
metric, call it the primary one. The other metric serves as an auxiliary to which the primary metric couples. Such a 
bimetric description of gravity has a long history, and, in particular, has been discussed recently in connection with 
some "dark substance" phenomena (see [16-19], and references therein). The novelty in BIMOND is in the special 
choice of interaction between the two metrics, which is particularly germane in the context of MOND, and which 
reproduces MOND phenomenology. 

In [15], I made preliminary comments on cosmological solutions within BIMOND. Classes of cosmological solutions 
within certain versions of BIMOND were also considered in [20]. But, there is still much left to be done on this 
problem. 

We are also still lacking any treatment of the important problem of the growth of small perturbations. Ideally, one 
should treat the problem in a fully relativistic framework, as done in [12-14], in the context of TeVeS-like theories. Here 
I take the lighter task of considering only density inhomogeneities of NR matter on scales smaller than cosmological. 
This limited-scope formalism, which governs the growth of perturbations under matter dominance, is still rather 
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illuminating, and shows a variety of phenomena potentially highly pertinent to structure formation in our universe. 

I limit myself to a universes that on average is completely symmetric in matter and twin matter (TM). Twin matter 
is the putative matter that couples to the "auxiliary" metric as (ordinary) matter couples to the "primary" one. At 
present, the existence of TM is not required phenomenologically (in particular, it is not, in fact cannot be, a sort of 
"dark matter," as assumed to explain galaxy dynamics); it is invoked mainly for aesthetic reasons: With TM we can 
obtain a simple and symmetric world picture that affords favorable cosmologies. In addition, as our present analysis 
shows, symmetric presence of matter and TM has very interesting consequences for structure growth. It leads to 
phenomena that are not only greatly different from those in GR (without dark matter), but also different from those 
in MOND without TM. 

Of course, the existence of TM would add crucial unknowns to our cosmological world view. Its properties, 
composition, self-interactions, amounts, etc., are not known (see [15] for more on this). To minimize these unknowns I 
assume here that, with one exception, the TM sector is identical to the matter sector in all respects, having the same 
properties (masses, etc.), the same interactions within the sector (electromagnetic, strong, etc.), the same interaction 
with its metric (geodesic motion, etc.), and the same average amounts of its different components present in the 
universe. As shown in [15] this symmetry leads to simple cosmological solutions, identical with those of GR, albeit 
with a naturally appearing cosmological constant (CC) of order (ao/c)^, where Uq is the MOND constant.^ The only 
difference between the two sectors I allow for, is in the amplitude and spectrum of the primordial inhomogeneities. 
These are putatively ascribed to quantum fluctuations in the early Universe. They could, thus, differ in the two 
sectors-at least in their exact distribution, if not in their statistical properties-even with all else being the same. 

Some preliminaries of the problem discussed here arc alluded to in [21], where I considered the problem of NR, 
small, isolated matter-TM systems. Here, this problem is extended to the evolution of perturbations on an otherwise 
homogeneous and matter-TM symmetric (MTMS), expanding universe. 

In section II, I briefly describe the aspects of BIMOND and of its NR limit that we need in this paper. In section III, 
I derive the equations controlling the cvohition of small departures from exact density equality in general NR systems. 
This is applied to cosmological small inhomogeneities in section IV. In section V, I describe succinctly how the density 
differences back react on the cosmological equations, introducing modifications that can act as cosmological "dark 
matter" and "dark energy" (beside the CC). Section VI is a discussion. 

II. BIMOND RECAPPED 

BIMOND is a class of relativistic theories of gravity that involve two metrics, g^^ and g^^. They are governed by 
an action of the form 

/=-Y^ jW'^R + ag''^R-2{ggfl^alM]d^x + I^{g,^,^l^i) + i^{g,^,i>i), (1) 

where R and R are the Ricci scalars of the two metrics as appear standardly in the Einstein-Hilbert action of GR, G 
is Newton's constant, and I use units where c = 1. The interaction term is a dimensionless, scalar function of the 
two metrics and their first derivatives. 

Bimetric theories have been much discussed before; the novelty and crux of BIMOND is in the choice of the 
interaction term. The key observation, in the MOND context, is that with two metrics and their first derivatives, we 
can construct tensors with the dimensions of acceleration: 

where F"^, F"^ are the Levi-Civita connections of the two metrics. Within the MOND paradigm, having the acceler- 
ation constant, Oq, at our disposal, we can construct from these dimensionless scalars to serve as variables on which 
M. depends. These scalars are obtained by contracting the dimensionless tensors 0,7 ^C"^.^ It befits us to concentrate 
on scalars that are quadratic in the C^^, which I do. In particular, the scalars constructed from the quadratic tensor 

= (3) 



^ So the successes of standard Big Bang cosniology, such as nucleosynthesis, are retained. The symmetry assumption entails, for example, 

having the same inflation, and the same mechanism for baryogcncsis in the two sectors. 
^ In the units where c = 1, ao has dimensions of length^ ^ . In standard units, ao is replaced everywhere by (.~^ = Oo/c^. Then in 

the NR limit only oo appears. In standard units, C^^ have dimensions of length~^; so the dimensionless tensors to be used are 
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have particular appeal; this tensor is symmetric under the interchange of the two metrics. The scalars I use are then 

T = 5'"'T^,, t = 5'"'T^,. (4) 

These are obtained from each other under metric interchange, and are used symmetrically in constructing MTMS 
actions for BIMOND.^ The action terms and Im are the matter actions for standard matter and for the putative 
TM (with their degrees of freedom denoted by ■0^ and ijji), whose existence is suggested by the double metric nature 
of the theory: Just as standard matter couples only to the standard metric g^^, TM couples only to the second metric 
There is no direct interaction between the two matter sectors. The dimcnsionless parameters fi and a permit 
us to use gravitational couplings in the two sectors, G' = GJ j5 and G' = G/a, that differ from each other and from 
G itself. In the rest of the paper I shall concentrate on symmetric versions of the theory; this means taking a = /3, 
and that the two metrics appear in a symmetric way in the interaction term. As explained in [15, 21], there are good 
reasons to take /? w 1 (see below). However, for the sake of generality I shall keep a general value of /3 in some of the 
discussion. Note that without the interaction term the theory consists of two uncoupled copies of GR. 
The resulting field equations are of the form 

pG,^ + S^^ = -87vGZ^, (5) 

pG^. + S,^ = -87rG%^, (6) 

where 

Gf,^ = Rf,^ - -jp^^R, G^^ = R^^ - -%^R (7) 

are the Einstein tensors, and 7^^,, are the energy-momentum tensors (EMTs) for the two sectors. 

The equations of motion of matter in the two sectors, in the presence of gravity, are the standard ones since they 
are still controlled by the same matter actions as in standard physics. 

The tensors S^^ and S^^ fall from the variation of the interaction term with respect to g^^ and respectively."^ 
The metrics appear in the interaction in three ways: (i) through their derivatives, in Cj^, (ii) contracting the C^^ in 
the quadratic scalars on which M depends, (iii) otherwise, such as in the volume density (gg)^^'^, or in scalars such 
as g/g or g'"'%^. As a result, S^,^ and S^,^ have the schematic form 

S^u = {Q{Crj.,. + N{CC}^. + alPg,., 

S^. = {Q{CVJ... + N{CC}^. + alPg,^, (8) 

with the three term types corresponding, respectively, to the above types of appearance of the metrics. Here, {C}^^ 
denotes tensors linear in Cj^,^ and {CG)^i, denotes tensors quadratic in them;^ they are the same in the two equations. 
Also, ":" denotes covariant derivation with respect to g^„. Q, N, P, etc. are dimcnsionless, and presumably of order 
unity. They depend on the two metrics through the quadratic scalars, and possibly through scalars not containing 
the metric derivatives, such as g/g or g^"g^^- Q and N are linear in derivatives of M. with respect to the quadratic 
scalars, e.g., T and T defined in Eq.(4), while P gets contributions from M itself and from its derivatives with respect 
to scalars containing the metrics alone. This classification of the tcirnis helps understand some interesting special 
cases and limits of the theory. In what follows I refer to these contributions as Q, N, P terms, respectively. 



^ The special appeal in the choice of the tensor T^i^, and the two scalars constructed from it, is twofold [15]; First, it leads to a theory 
with a particularly simple NR limit. Second, this choice is the best equivalent to the choice of the Ricci tensor, Rpiv, and the Ricci 
scalar, R = g''" R^v, in the Einstein-Hilbert Lagrangian of GR. In GR we have to use the full curvature tensor, which contains the 
second derivatives of the metric; this is both necessary (to get a scalar), and harmless (because the second derivatives in the Einstein- 
Hilbert Lagrangian do not affect the field equations). In iaict, it is well known that we can repla<;e ii in the GR Lagrangian by 
9^" " '^/Hi-'^At)' ^^''^tly analogous to our T, T scalars. In our case we do not add the second derivatives; this would be both 
harmful (as it would lead to a higher-derivative theory), and unnecessary (as the first-derivative terms are tensor in themselves). 

* They satisfy, together, one set of four Bianchi-like differential identities, stemming from the invariance of the interaction term to general 
coordinate transformations (details in [15]). 

^ Examples of such tensors are C^„, S^C°^, gi^i^g°^C^i^, etc. 

® Examples of such tensors are C^xC^^, C^^C^^, T^i,, etc. 
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Prom the symmetry of the problem we have 

Q{9^..,%u) = -Q{%.,9^..), N{g^,,g^,) = N{g^,,g^,), P{g^.,%S) = P{%.,9^S)- (9) 

BIMOND versions with (3 = 1 have a simple GR limit: Assume that in the limit Oq 0, i.e., when the quadratic 
scalar arguments of M. go to infinity, M. becomes independent of them. Then, the Q and A'' terms vanish. The P terms 
then give CC-type terms that might still couple the two metrics, but which are rather small locally. It is also possible 
to choose Jvl such that the coupling disappears altogether in this limit, so the P terms in S^^ collect to a^P(oo)g,^^, 
[15]. The theory then decouples into two copies of GR, possibly with a cosmological constant, if A^(oo) ^ 0. I have 
not yet been able to ascertain if there is an acceptable GR limit for theories with (3 ^ l7 This may be one reason to 
prefer theories with /3 = 1. 

When the two metrics are equal, we have from Eq.(9) Q = —Q = — Qo, N = N = No, P = P = Pq, where, 
furthermore, all these functions of the two metrics become constants. This is because we cannot construct (non- 
constant) scalars from a single metric and its first derivatives, so all the scalar arguments of these functions must 
become constants when = g^^. Clearly, Q terms vanish linearly in the metric difference, and terms of type N 
vanish quadratically. We are then left only with P terms, which are of a cosmological constant type with A oc a^Po. 
From this follows that, in systems with identical matter and TM distributions, and the same initial and boundary 
conditions, the metrics in the two sectors are the same, and are a solution of the standard Einstein equations with a 
cosmological constant (CC) A.^'^ 

A. Cosmological background 

Preliminary considerations of cosmology in BIMOND have been described in [21] and in [20]. As mentioned in 

the Introduction, here I assume all along the special, but attractive, case of an MTMS cosmology. This means that 
matter sources in cosmology are assumed to be identical in the two sectors, apart from differences in the initial seed 
inhomogeneities. If we neglect these differences to lowest order, the MTMS cosmology in BIMOND is identical to the 
Friedmann-Robertson- Walker (PRW) cosmology in GR, with a CC, A cx a^Po, with G' = G/fH instead of G, and with 
the standard matter content. All the standard results derived for the latter (Big Bang, Nucleosynthesis, etc.) apply 
here as well (with G' as gravitational constant). 

To a higher order, the differences in the inhomogeneities in the two sectors engender departure of the background 
cosmology from PRW: the cosmological metric is no longer a solution of the Einstein equations. This is discussed 
briefly in section V. To some approximation the cosmology could still be considered an FRW one, but with non- 
standard, "phantom" contributions to the cosmological EMT in BIMOND. This sort of feedback goes beyond the 
feedback due to departure from smoothness, which has been discussed extensively respecting GR cosmology. Here it 
is specifically due to differences between the inhomogeneities in the two sectors. As long as the statistical properties 
of the inhomogeneities in the two sectors are the same, the cosmological background metrics are still the same, but 
they satisfy the BIMOND equations, instead of the Einstein equations. 

The observed "^He abundance was used in [22] to derive a bound of \G'/G — 1| ^ 0.13, where G' is the value that 
applies at the formation time. This would imply, in the present context, |/S — 1| ^ 0.13, lending further motivation 
for taking ^ = 1. 

B. The nonrelativistic limit 

Consider first the NR limit of BIMOND as applied to a system that is small on cosmological scales (discussed in 
detail in [15]). Since we assume a symmetric cosmology, with the two metrics equal, there is a coordinate frame in which 
locally the two metrics are nearly Minkowskian. Because of local inequality of matter and TM mass distributions, 
the departure from Minkowski is not the same for the two metrics, which can thus be written as 

9i.. = r)i,^-2(j)5^^ + hu,^, g,,^ = f]^^ - 2(j)6^^ + hf,^, (10) 



The NR limit of the theory does have a good Newtonian limit even for ^ ^ 1. 
* For example, a double Schwarzschild solution is a vacuum solution of BIMOND, which corresponds to a black hole of equal matter and 

TM masses. Similarly, standard GR gravitational waves are also vacuum solutions of BIMOND. 
^ Pictorially, in a double membrane picture the Einstein-Hilbert actions for the two metrics may be viewed, as usual, as the elastic energies 
of the membranes— the cost of distorting them out of the flat. Then, the CC may be viewed as the contact energy density, remaining 
when gft„ = g^i,, and the rest of the interaction is the "cost of separating" the membranes. 
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with the potentials cj), cj), h^^, h^i, treated to lowest order in the weak- field approximation.-'^'^ 

Note, importantly, that the weak- field approximation assumes that </>, <C 1 (in our units where c = 1). This, 
however, does not imply that whenever <p or (j) appear they raise the order in the perturbation: Unlike the case of 
GR, here we have an additional dimensioned constant, a^, and quantities such as V(j)/ao or (V<^ — V(f))/aQ are not 
assumed small, and do not raise the order of terms in which they appear. In this light, as they appear in {C}^^ 
and {CC}^i, are first order in the small parameters; however, in their appearance in the arguments of Q, N, P, etc. 
cin^Ca should not be taken as small. -'^^ 

The N terms are second order and can be neglected. The P terms contribute to lowest order as a homogeneous 
CC, whose value we know is observationally negligible in local NR physics. As to the Q terms, since the {C}^^ are 
already first order we can equate the two metrics everywhere else with the background Minkowski metric (so, for 
example, covariant derivatives become normal derivatives). But note that this approximation applies only to the 
metrics themselves, not to their derivatives, which appear only in a~^C^^, which are not assumed small. 

It can be seen^^ that to the lowest order Q = —Q, and so = —5^^. Thus, the sum and difference of the two 
field equations are, in this approximation, 

PiG^^+G^^) = -8nG{%^+%^), p{G^^-G^^) + 2S^^ = -8nG{%^-%^). (11) 

The sum of the Einstein tensors is linear in the sums of the potentials, and their difference, as well as 5^^, depends 
only on the differences of the potentials. 

In [15] I treated the NR limit where one further assumes that motions are slow, so that time derivatives are negligible 
compared with spatial ones. I then showed that with the choice of the quadratic scalar arguments given in Eq.(4), 
which I shall assume, the field equations imply, in a certain choice of gauge: h^^, = /i^i, = 0. The two metrics thus 
have the form of the first-order metric in GR, but the NR potentials, </> and (j), are determined not from the Poisson 
equation but from field equations derived from the action 

£ = -^{/3(W)' + /3(V^)^ - alM[{Vct> - ^4>r/am + p{^v' - <^) + p{h^ - (12) 

This follows from the fact that to lowest order we have T = T oc (V(/) — V(/))^, and so the interaction term takes the 
form a^Mliycj) — Vcf))^ /a^]. (I use here the same symbol, M., for a function of a single variable: the reduction of 
the relativistic interaction function, also denoted M, which depends on more variables.) This action also governs the 
particle equations of motions in the two sectors, which take the standard form: a = — V0 and a = — V(/). Because the 
metric has the GR form in terms of a single potential: = ?7^„ — 20(5^„, light bending by slowly moving masses is 
expressed in terms of the potential in the same way as in GR: photons "see" the same potential as massive particles. 
The field equations for the potentials are 

= 47rG/3-V + V • [(V^ - ^4>)M'] = 4nG{p + pp), 



A4> = A-KGp-^p - V • [(V0 - W)-A4'] = 47rG(p + Pp), (13) 

where pp and pp play the role of "phantom matter" (PM) densities for the two sectors. 
As detailed in [21], it is convenient to work with the two potentials ^ and such that 

<p = C$ + l ^ = C^-4>, (14) 

which satisfy 

A0 = 47rG(pH-p), V-{M(|WI/ao)W}=47rG(p-p), (15) 



The separation into (p Etnd h/xv is a mere convenience: One defines (j> = —(1 -\- poo)/2, and hoo = 0. 

Oo may be viewed as a proxy for a scale length i = c^/oq. Terms such as V</>/ao appear, in fact, as ^V(</>/c^). The NR limit in the 
MOND context has to be understood as c — >■ oo, t — >■ oo, with c^/^ fixed. 

Following from Eq.(9), noting that in this approximation Q is symmetric to the interchange of the metrics, as its argument is quadratic 
in their difference. 

The general form of the scalar was also discussed in [15]. 
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with the boundary conditions at infinity (j) ^ 0, while (j) ^ ln{r) when M =^ M, and — > when M = M; 
M and M are the total masses in matter and TM respectively. Here, ( = (2/3)~^, and fi{x) = 2/3 — 4A^'(4x^).-^'^ 
MOND and Newtonian behaviors for isolated, pure-matter systems, such as galaxies, dictate the asymptotic behaviors 
jl{x <C 1) « X, jl{x ^ 1) sa (1 — Thus (f) = (f)'^ +({)" is simply the sum of the matter and TM Newtonian potentials, 
and ^ is a MOND potential satisfying the nonlinear Poisson equation proposed by [23] with p — p as a source (and a 
modified value of at large arguments). 



III. FLUCTUATIONS IN AN ISOLATED MTMS SYSTEM 



Consider first the development of weak fluctuations in a cosmologically small, isolated, NR system. The two matter 
components are treated as fluids. The time dependent problem is governed by the NR continuity equations in the 
two sectors 

+ V • (pv) =0, p, + V ■ (pv) = 0, (16) 

and the respective Euler equations 

Vt + (v • V)v = -V0 - p"'Vp, Vt + (v • V)v = -V(^ - p"'Vp, (17) 

where p, p are the pressures in the two sectors. 

In a system made predominantly of matter-as we have, purportedly, in the solar system, galaxies, etc. -these 
equations lead to standard MOND dynamics as they have been applied extensively in the past. In this case, the 
MOND acceleration is generically larger than the Newtonian one, and when the acceleration is smaller than Cq, 
dynamics is controlled by the MOND limit. 

However, as discussed in [21], things are rather different in systems where the distribution of matter and TM are 
close to each other (i.e., |p — p| -C p). For example, when p = p the dynamics is always Newtonian, no matter how 
small the acceleration is.^^ 

Consider then weak perturbations on an arbitrary, symmetric background flow with p = p = P(,(r, t), v = v = 
v^{r,t). Then, from Eq.(15), ^ = 0, and = = /3"\6f , where 0" is the Newtonian potential of p,,: we have a 
double Newtonian system with an effective Ge = G / (3}^ If we seed the system with fluctuations that are the same 
in the two sectors, they will remain so, with = at all times, and the development remaining quasi-Newtonian in 
both sectors. 

Taking the lowest order in the fluctuations we have 

5p, + V • {5pv,) + V • (p,(5v) = 0, 5p, + V • ((5pv,) + V • {p,5w) = 0, (18) 
and the Euler equations 

6v, + {6v ■ V)v, + (v, • V)(5v = -CV54> - V0 - p^'VSp + 5ppl^Vp„ 

5w, + {5w ■ V)v, + (v, • V)(5v = -CV50 + V0 - P,"' + Spp^'^Vp,, (19) 

where pt is the pressure field in the unperturbed system (assumed to be the same in the two sectors, as we assume 
the same equation of state for the two). The potential fluctuations are obtained from: 

A64> = 4'KG{6p + Sp), 

V • (/i I V0/ao I V0) = ATrG{Sp - 6p) . (20) 
If \5p\ ~ \5p\ ~ \5p — 5p\, the sources for these equations are of similar magnitudes; but this need not be so. 



The factor 4 in the argument accounts for the fact that we use the variable (f) = ((j> ~ (f)) /2. 

This is the NR expression of the reduction of BIMOND to GR for completely symmetric systems. 

This case is not to be confused with the Newtonian limit in a system made of pure matter or pure TM, in which case G itself plays the 
role of the Newton constant. 
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It is useful to work with the sums and differences: 

u = + Sir, u = 6v — Sir, q = 6p + 6p, q = 6p — 6p. (21) 

If we can assume that in both sectors p can be taken a function of p; e.g., for a polytropic gas, or when the process 
is isentropic, we write 5p = (dp/ dp){pb)5p = vl{pb)5p and get: 

+ V • {qw,) + V • (p,u) = 0, (22) 

u, + (u • V)v, + (v, • V)u = -2CV50 - p;'V{vlq) + qp-,'vyp,, (23) 

A(5(^ = 47rGg, (24) 

+ V • (gv,) + V • (p,u) = 0, (25) 

u, + (u ■ V)v, + (v, • V)u = -2V^ - p;'W{vlq) + qp^'vyp,, (26) 

V • (/i|V0/aolV<^) = AirGq. (27) 

The barred and unbarred systems decouple, with the unbarred fluctuations developing according to standard New- 
tonian theory, while the barred fluctuations develope according to MOND. If we start with q and q fluctuations of 
a similar magnitude, and the barred fluctuations are in the MOND regime, q grows faster than q. This leads in the 
limit to Sp « —Sp, namely segregation of the two sectors. This echoes the fact that matter and TM repel each other 
in the MOND regime [21]. 

IV. COSMOLOGICAL FLUCTUATION GROWTH 

The use in the cosmological context of the above NR equations, derived for an isolated system on a double Minkowski 
background, can be justified for weak perturbations in NR matter, of scale length much smaller than cosmological 
scales: Expand the BIMOND field equations around their cosmological solution of the FRW equations with a CC and 
9iiv = g^v The EMTs are expanded around their value for the homogeneous background with only the 00 components 
contributing. For scales much smaller than the (space-time) curvature radius the background metrics can be taken as 
nearly Minkowski. Then, our formalism above applies. 

In this application, pt{t) is the ambient, homogeneous cosmological density of the specific matter component we 
concentrate on, and pt. is its effective pressure. 

I shall take in what follows (3 = 1. 

As is standard, we express the space derivatives with respect to comoving coordinates R instead of the proper- 
distance coordinates r — a(t)R, which appear in the above equations. We also put v^ = H{t)r, V ■ — 3H, where 
H = a/a is the Hubble parameter. Also, write the equations in terms of partial time derivatives at a fixed comoving 
position: for the field S, S = St + (v^ • V)^. The continuity equation for the background NR matter amounts to 

oc a-"", so, p,, = dpt/dt = -3Hp^. 

Define r] = 6p/pi,, e = t] + i] = q/ p^, e = r] — r) = q/ Pb- For these we have: 

e-|-a-'VR-u = 0, (28) 

u + (d/a)u = -a-^V^5(i> - vla'^V^e, (29) 

AaJ^ = A-KGa'p^e, (30) 

e-|-a-'VR-u = 0, (31) 

u+ (d/a)u= -2a-'VR^ V^e, (32) 
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VR(/t|VR^/aao|VR(^) = 47rGaV,e. (33) 

The unbarred system is identical to that in standard gravity: the perturbation sum develops fully in a Newtonian 
manner. Combining the time derivative of the continuity equation with the R divergence of the Euler equation we 
get as usual 

e + 2(d/a)e = AnGp^e + v^a'^ l^^e, (34) 

and for the space Fourier components defined by e = / cPkey_{t)e^^ '' , we have the standard result 

e'k + 2(d/a)ek = (47rGp, - vlk^/a')e^. (35) 

The barred system, which is decoupled from the unbarred one, describes the evolution of the perturbation difference; 
it behaves as a single fluid governed by MOND, with the strength of the governing potential, cj), doubled [as expressed 
by the factor 2 in front of V0 in Eq.(32)], and with a Newtonian limit for which the effective Newton constant is G.^^ 
We can still eliminate the velocity u to get 

e + 2(d/o)e - 2o-' Ar^ - ArC = 0, (36) 

but here Ar^ is not simply expressible in terms of e. The problem is nonlinear even for the smallest of perturbations, 
and perturbations on different scales can affect each other's growth (see below). 

_ Given some initial e(R, to) and u(R, io), we can calculate e(R, io) from the continuity equation, we then calculate 
(j) from the potential equation, then calculate e from Eq.(36), and so propagate the system in time. The velocities are 
propagated in time in parallel, using the Euler equation. 



A. The deep-MOND regime 

With matter and TM present, there are two senses in which we can understand the term "a deep-MOND system": 
The weaker sense applies to systems where [Vet)] <C cia so that is determined by deep-MOND physics (with its 
peculiar symmetries, etc.). If we have an inhomogeneity of proper length scale A that can be considered in isolation, 
i.e., is not subject to an external- field effect (EFE) from a larger-scale inhomogeneity of which it is a part (see below), 
this amounts to 

\5p - 5p\G\ < ao- (37) 

In the stronger sense of the term, one requires, also, that |V(^| <C |V^| so that the full potentials <f) and (j) are 
determined by deep-MOND physics. This requires further 

\5p + 8p\G\ < {a^\8p - 5p\GXf'^. (38) 

In past applications of MOND, where TM is absent, the second requirement follows automatically from the first. This 
is also the case in the present context if \5p — 5p\ ^ ■max{\5p\, \5p\) (for example, if the fluctuations in the two sectors 
are uncorrelated) , but not if \5p — 5p\ <C max{\5p\, \5p\). 

Because here we have full decoupling of the sum and difference systems, with the latter being governed only by </>, I 
shall use the term in the weaker sense, which is enough to ensure that the barred system is governed by deep-MOND 
physics. 

In the deep-MOND limit, Eq.(33) is well approximated by 

Vr(|Vr^|Vr^) = 47rGpoaoe, (39) 
where po = a^Pb- This equation is conformally invariant [21, 24]. 



Since p-{oo) = 2, is half the Newtonian potential sourced by pbe in this limit; the 1/2 factor canceled by above mentioned factor 2. 
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B. Coupling of inhomogeneities on different scales 

The impUcations of Eq.(33) for the dynamics in massive systems have been discussed extensively (e.g., in [23, 25]). 
One clear outcome of these studies concerns the way in which the dynamics of a larger, mother system coiiple to those 
in a smaller subsystem. When the two do not differ much in size and/or mass, the interaction between them, which 
results from the nonlinearity of the theory, is not easy to describe. However, when the subsystem is much smaller, 
and much less massive, than the mother system (such as a for star in a galaxy, or a galaxy in a galaxy cluster), 
the interaction is simply described: The subsystem moves within the mother system as if it where a test particle, 
irrespective of its own intrinsic accelerations. In other words, the dynamics of the mother system is oblivious to the 
internal dynamics in its small subsystems, which can be treated as test particles moving in their combined mean field. 
On the other hand, the dynamics within the subsystem may be greatly affected by the acceleration, g^^, with which 
it falls within the mother system, provided this acceleration dominates the internal ones, and is in the; MOND regime 
(ffex ^ do)- The intrinsic dynamics is then quasi-Newtonian with an effective gravitational constant G' ~ Guo/g^^, 
and some mild anisotropy induced by the external field. This is known as the MOND external-field-effect (EFE). 

In the present context, the mass source in the equation for the potential, must exhibit, at any time, a continuum 
spectrum of fluctuations, and it is difficult to describe the mutual couplings of the different scales. We expect, however, 
that the evolution of a fluctuation in e on a given scale, will be affected mainly by larger-scale perturbations within 
which it is nested, in the spirit of the EFE. 

A difference fluctuation of magnitude e and (proper) scale length A, would have a characteristic intrinsic acceleration, 
when in isolation, ^ {agGpteXy^^. Thus, if, e.g., e oc A"°, for a > 1 the typical acceleration decreases with increasing 
scale, and we do not expect a pronounced EFE. For a < 1, the largest scales determine the dynamics within smaller 
ones, as described roughly by the EFE. For a = 1 the accelerations on all scales are of the same order. 



C. Pressure vs gravity 

For the unbarred, fluctuation-sum system, the standard (Jeans) criterion tells us when pressure becomes unim- 
portant in impeding gravitational growth. This applies [e.g., per Eq.(35)], when p^GX^ ^ v^, or, equivalently, when 
Mb ;S> G~''''^vlpy,^'^ , where M,, = pi,\^ is the background mass (in either sector) in a cube of side A. 

For the barred system, the Jeans criterion for gravity dominance is very different. We saw that because of the 
nonlinearity of this system we cannot consider separately the growth of a single perturbation with a given amplitude 
and scale length, as all modes are, in principle, coupled. I discuss two end cases. First, consider a perturbation of 
magnitude e and (proper) scale A that is not affected by larger-scale overdensities (through an analog of the EFE). 
From Eq.(32), or Eq.(36), the Jeans criterion is still |V^| > I'JjVel. Prom Eq.(33) we now have in the MOND limit 
[V^l^/Aflo ^ Gph\e\. Put together, these give the condition for gravity dominance 

M, > {Gao)-X\e\ (40) 

(where gradients are substituted by division by A; c.g.,|Ve| ^ jej/^)- 

Note, in particular, the appearance of |e| in the criterion, in addition to the expected characteristic, deep-MOND 
M oc V* relation. This occurs because the pressure increment is linear in e, while the potential increment 4> is 
proportional to |e|^'^. This |e| factor makes the criterion much easier to satisfy for weak perturbations, -'^^ 

If the particular perturbation at hand cannot be considered in isolation, but is subject to an EFE characterized by 
5ex, tlie Jeans criterion takes the form p^G{ao/ gex)X? '> v^, or 

M, » G-'''{gJa,rXp;'''- (41) 

As a concrete example, consider a baryonic, difference perturbation, of magnitude e and scale A, that may be 
considered unaffected by an EFE. Taking the ambient temperature appropriate for the matter-dominance era T w 
3 X 10^(a/I0"^)"^K, and ambient density « 4 x 10'^^(a/10"^)"^gr cm~^, the condition (37) for the perturbation 
being in the deep-MOND regime can be written 

Mb/Me «2x 10^ (^)'|e-|-. (42) 



In [26] I gave the Jeans criterion for a MOND universe as Mb ^ (Gao)^^«f . This was based on a pure matter world- picture, where 
it was assumed, tentatively, that the source for the MOND potential includes the background density, and so the potential increment 
is linear in the density increment. Here, with an MTMS background, the background MOND potential vanishes, so the perturbation 
alone sources 0, and so |V^| oc 
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The condition (40) for gravity dominance over pressure is 

M,/A4»3xlO=(^)"|6-|. (43) 
In comparison, the Jeans criterion relevant for the sum of the perturbations is 

M./Me»10«(^)"". (44) 
So, there is a wide range of scales for which pressure stymies the growth of e but not of e. 

D. Later stages 

The analysis above highlights the way in which small-scale, NR density inhomogeneities grow when they are still 
small compared with the ambient density. We have learned much about the mechanisms that are at play and how 

they greatly differ from those underlying perturbation growth in standard dynamics. For example, we noted the 
shepherding of matter overdensities by flanking TM overdensities (and vice- versa) : a matter overdensity grows faster 
when it sits in a TM depression. A related effect is the segregation of matter and TM, as evinced by the tendency of 
difference fluctuations to grow faster than the sum 

But, beyond these qualitative observations, a numerical attack on the problem is needed to follow the growth of 
even weak perturbations in detail, because of the inherent MOND nonlinearity. This is also, clearly, the case when we 
deal with structiirc buildup to its presently observed stage, where the perturbations arc not small. As mentioned in 
the introduction there have been several A^-body studies of structure formation in MOND, and these could be easily 
adapted to include TM. 

Without numerical studies, it is difficult to tell exactly how, and to what extent, these trends continue to shape 
structure in the later stages. Some questions that naturally arise are: "Up to what scales are matter and TM effectively 
segregated today?" In [21], I conjectured that this happens on the supcrcluster-void scales; namely, that matter and 
TM have formed interlacing cosmic webs in which TM concentrations sit in the observed matter voids, and vice- versa, 
and matter and TM filaments avoid each other. Another important question is: "How rare are TM objects, such 
as galaxies, in matter territory, and what observable effects can they have?" If the emptiness of matter voids is an 
indication, we do not expect TM galaxies within large matter concentrations. Other questions are: "Does TM help 
empty matter voids more efficiently than in the CDM paradigm?" (For possible indications for unexpectedly large 
and empty voids sec [27-33].) "Arc higher very-large-scale velocities produced in the present picture than in the 
framework of CDM?" (See [34 36] for possible indications of such unexpectedly high velocities.) "Are objects with 
high peculiar velocities more common in our world picture than in the CDM paradigm?" (See [37] for a discussion of 
the issue in connection with the bullet cluster, and [38] for a recent assessment, and for references to earlier work.) 

We can get some guidance in pondering over these questions by considering the interaction between well separated 
matter and TM bodies. This question has been considered in detail in [21], which deals, in a sense, with the opposite 
of the situation we treat here of small departures from well mixed configurations. 

The basic result found there is that matter and TM bodies repel each other in the MOND regime, and (for /3 = 1) 
they do not interact in the high-acceleration regime. For example, using the conformal invariance of the deep-MOND 
limit of the theory, I calculated the force between two masses Mi and in the deep-MOND regime to be (defined 
as negative for attraction) 

F = --i^£^l^[(Mi ± M^f'^ - Ml" - Ml'\ (45) 
3 V 

where the plus sign applies when the two masses are in the same sector (in which case F is always attractive), and 
the minus sign when they are different (in which case the force is always repulsive); r is the separation. We see that 
for equal masses the repulsive force between matter and TM masses M is (4/3)(aoG)^^^M^^^/r, which is ~ 2.4 times 
larger than the attractive force between two masses M in the same sector. This can give us an idea of how far from 
each other separated bodies of matter and TM can separate over the Hubble time.^^ 

The characteristic separation speed for matter and TM masses M is ~ 2.3(aoGM)'/'' ~ 260(M/10'°MQ)'/*km s"^ 
In a Hubble time his produces a separation of several megaparsecs for galactic mass objects and a few tens of 



This is only a rough estimates, since such bodies are hardly ever isolated on the relevant scales. 



11 



megaparsecs for galaxy-cluster mass objects. A test mass in one sector in the (MOND) field of a mass M of the other 
sector has typical separation speeds about half of that given above: a massive TM cluster can shoot out a matter 
bullet with a speed of several thousand km s"-^.^*^ 

Another potentially useful result for a late-stage system made of separate matter and TM bodies is the general, 
'virial' relation that holds exactly in the deep-MOND regime: For a system made of matter and TM pointlike masses 
rrii and m^, respectively, at positions (fc running over all masses), we have 

k i j i 3 

where is the force on particle k. (The choice of origin is immaterial as ^ = 0.) For a system with equal total 

masses, the first term drops, which gives an expression that is always positive; this means that on average the forces 
are pointing outward, bespeaking a mean repulsion. For example, in a system made of N equal matter masses, m, 
and N equal TM masses, m, we get 

Er.-F. = -(aoG)^/^7Vm=*/^ (47) 

fc = l 

Thus, for the accelerations we have 

{n.a,) = ^{aoGmy'\ (48) 

The acceleration thus points outward on average, implying global self repulsion. The magnitude of the typical 
acceleration is ^ {aoGniY^'^ / R, where R is the characteristic system size. This is smaller by a factor ^ N^^'^ than the 
typical MOND acceleration in a pure-matter system of the same mass. Note, however, that this mean acceleration is 
only relevant if the matter and TM masses are very well mixed. If, in fact, they form separate matter and TM clumps 
within the system, as they would tend to do, these clumps would separate with higher accelerations. 



V. PHANTOM MATTER AND ITS BACK REACTION ON COSMOLOGY 



There has been much discussion recently regarding back reaction of inhomogeneities on cosmological evolution in 

GR. The issue concerns the extent to which we err by applying FRW theory (which assumes homogeneous energy 
density) to the mean cosmological energy density, instead of calculating cosmological evolution correctly with inho- 
mogeneously distributed energy, and then averaging the outcome. (For recent discussion of this issue, with references 
to earlier work, see, e.g., [39, 40].) The exact importance of this effect is still moot. 

Obviously, such effects exist in BIMOND as well, but they are not the subject here. In a MTMS universe governed 
by BIMOND there arc added back-reaction effects of inhomogeneities, which, as we saw, induce the appearance of 
phantom matter components that should have influence on cosmology. 

In a universe with different matter and TM inhomogeneities, the S^^ and terms in the BIMOND field equations 
(5) (6) do not vanish; so, what is their effect on the cosmological background? Strictly speaking, we then have to solve 
the full BIMOND cosmological equations. However, to some approximation, we may treat the resulting departure 
from GR perturbatively, by viewing the effects of these terms as perturbations on the MTMS, FRW cosmology: In 
this approximation, the S^^ and terms are viewed as the EMTs of MOND phantom matter that can effectively 
be added to and in sourcing cosmic gravity. We saw that to zeroth order in the metric difference, all that 
remains from these terms is a CC term, which has to be reckoned with in zeroth-order cosmology. This cosmology 
is then used to calculate the development of perturbations in the two sectors, starting from some seed perturbations. 
The next stage in such an approximation scheme would be to use the deduced local departures of the two metrics 
from the FRW metric, used in the lower order approximation, to calculate local contributions to S',,,^ and S^,^, and 
put them back as sources in the equations for cosmology. This program clearly deserves a thorough treatment, which, 
however, is beyond our scope here. In particular, the results of such an analysis is expected to depend on the details of 



If the bullet test mass starts at rest a distance Tq from a TM mass M, its speed at a distance r is, assuming the system is isolated, 

v{r) = (A,aaGMY'Hn^''^{r/ro) ~ 1.6 X W^{M /W^* MqY' Hn^''^{r /ro)^m. s"!. 

Such a system docs not satisfy a kinematic virial relation: since = mkVk, we find that (i(X)ic ^ki^k ■ 'Vk)/dt > 0; so no steady state 

can be reached. 
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the specific BIA/[OND version at hand, and there is a large variety of versions to explore. Here I only describe briefly 
some general properties of the different term types defined in Eq.(8), as they would be fed back into the cosmological 
equations. 

The zeroth order of the P term (in the matter sector), a^Poi?^^ gives a genuine CC term [if Pq 7^ 0]. When matter 
and TM separate, P becomes space-time dependent, and in this regard behaves as a "dark energy" field. In the weak- 
field approximation, we can substitute in the variables on which P depends w g^^, but T T oc (V^ — V^)^. 
The contribution of the P term is then 

« aoP{z)g,., z = {^cl>- V^f/al. (49) 

Prom Eq.(9) we see that P = P in this approximation. The value of P in this expression varies in space-time between 

Po and P(oc). (As explained in [21], considerations of the GR limit of the theory, for Oq 0, dictate that P(oo) is 
not infinite.) Today we have in much of the volume of the universe z < 10"^; so in much of the volume P ~ Pq, but in 
galaxies or in galaxy clusters (within about a megaparsec of the center), z>l, and near stars P w P(oo). Since the 
values of P around matter concentrations are the same as its values around similar TM concentrations (interchanging 
matter and TM amounts to interchanging </> and (j) to which P is invariant), the large-scale average of P gives a time 
dependent effective CC, A(t) oc (P),p„ee 7^ -Po- 

The next order contribution, which appears with the separation of matter from TM, comes from the lowest order 
in the Q terms, whose contribution to the cosmological density "phantom" EMT is 

r::l' = {inG)-\Q{CYJ,,, (50) 

giving in the weak-field approximation 

T;^' « -tr^ « (87rG)-(Q{C}^J,,. (51) 

In the weak-field-slow-motions approximation, where time derivatives are neglected, only the 00 component con- 
tributes and gives, simply, the totality of all the PM gotten from Eq.(13), namely a NR PM density 

Pp = (47rG)-'V • [(V^ - V^)M'{z)]. (52) 

In a MTMS universe this density averages to zero over large cosmological volumes. This can be seen in two (inde- 
pendent) ways. First, it follows from the fact that under the interchange of matter and TM, 7^^ O 7^„, we have 
g^u^ ^ In the NR limit this means that p p leads to (/> -H- </>, which, in turns gives Pp — >■ —pp. The phantom 
density around a matter concentration is equal and opposite in sign to that around a similar TM concentration. In 
our world picture the large-scale averages of the universe are invariant to matter-TM interchange, hence the average 
phantom density, which changes sign, must vanish. This result holds also for a finite system with matter-TM symme- 
try. For example, in a system where matter-TM interchange can be achieved by a volume-preserving transformation 
(such as translation, refiection, rotation) the total phantom mass in a symmetric volume vanishes. 

We can also argue differently, noting that the phantom density is a divergence of a vector (V^ — V(/))A^'(2), whose 
absolute value is bounded (outside of point masses). As a result, its volume average vanishes in the limit of a large 
volume. This property of the phantom density follows from the form of the relativistic Q term contribution: We see 
from Eq.(50) that the trace of the phantom EMT is a covariant divergence of a vector oc <3i?'"'{C'}^^- Its integral 
over a space-time volume is hence a surface integral. 

Even if the large-scale average of the NR phantom density averages to zero it could still affect cosmology due to its 
inhomogencity. In addition, there may be higher order contributions that do not average to zero. 

The Q terms may also have contributions of the form (A^'ff°'^C'^^),A(/^i,, or [M' g'^^C^-^.^fig^^, which are of the form 
V^-xg^,^; i.e., they contribute in cosmology as dark energy with oscillating, vanishing-average density. 

The PM around matter concentrations can be (and is) observed, and is interpreted by most as "dark matter." As 
discussed in [21], the PM around TM concentrations may be observed with weak lensing of matter photons (it does 
not produce strong lensing), on which it acts as a repelling mass. 

For the N terms, the contribution of lowest order in the departure from metric equality is the second order. In the 
NR limit their characteristic contribution to the phantom EMT is of the order of {^-kG)^^ zN [z)(aal cf [z defined in 
Eq.(49)]. Using the well known proximity Qq « cHo/2n, we can write this contribution as ~ 10~^zN{z)p^, where p^ is 
the critical density today. We saw above that in most of the cosmological volume today z < 10~^. Also we have that 
A''(0) is of order imity (from considerations of the MOND limit, see [15, 21]), and 7V(oo) = (from considerations 
of the GR limit of BIMOND). N is related to the extrapolating function of MOND and is constrained to remain of 
order unity, or below, for all arguments. So the contribution of the N terms to the density can be at most of order p^ 
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(in systems where |V(/)| ~ Oo), but some orders of magnitude smaller than in most of the cosmic volume. It would 
be interesting to see what their tensorial form is, and to what extent they prodiice observable lensing effects. 

Note, finally, that inasmuch as the statistical properties of matter in the two sectors are the same, so are the 
feedback effects discussed above. Namely, the effects of inhomogeneities, and the statistical properties of the PM 
(their cosmological mean densities, power spectra, etc.) are the same in the two sectors. The two mean background 
metrics thus remain the same. 



VI. DISCUSSION 

I have considered the dynamics of small, NR perturbations in an MTMS universe. An intriguing aspect of this 
world picture is that MOND dynamics, cosmological or local, are completely absent in a perfectly symmetric universe, 
apart from a possible appearance of a CC of order (ao/c)^. All other departure from standard physics encapsulated in 
MOND-as it appears in cosmology, in structure formation, and in bound systems, such as the solar system or galactic 
systems-is entirely an effect of "separation" between matter and TM, or separation between the two metrics. In a 
double-membrane picture of BIMOND, MOND effects are due to local separation of the membranes. 

Twin matter plays two roles in our narrative: In the first, it provide a counteracting medium to matter; similar to 
the role of a neutralizing positive ion background to an electron gas. As a result of this presence, we could assume the 
very appealing symmetric cosmology, identical to that of GR with a CC. In addition, the presence of matter and TM 
in equal average quantities leads to a well-defined MOND dynamics, in which the MOND potcintial is sourced only 
by departures from uniformity, not by the bulk. To fulfill this role, TM only has to be present in the same amount as 
matter; it does not necessarily have to clump. 

In its second role, TM plays an equivalent part to that of matter; so it can be likened to a positron component in a 
globally neutral electron-positron plasma (with mutual repulsion, and no annihilation). In this role it also takes part 
in shaping the large-scale structure of matter. 

There are several qualitative features that make the growth of fluctuations in the proposed world picture different 
from the standard one. Among them: (i) Two initial fluctuation spectra, not one, determine the large-scale structure 
at later times. For the NR problem discussed here, we can think of these as the fluctuations in matter alone, and 
of those in the matter- TM density difference. This is important to appreciate, for example, because we have strong 
constraints on matter fluctuations at the time of matter-radiation decoupling, from the CMB anisotropies; but, we 
have none on the difference fluctuations e, at that time. So, numerical simulations have to explore all possible options, 
(ii) The NR fluctuations in the sum and difference of the densities form two decoupled systems controlled by different 
types of gravitational physics. The sum system is underlaid by standard Newtonian gravity, while the difference 
system is underlaid by MOND gravity, which is stronger for the same source strength, (iii) In the difference system, 
the competition between gravity and pressure is strongly tilted in favor of gravity, with the Jeans criterion being 
much more lenient than in the sum system, (iv) As a result of (ii) and (iii), the difference fluctuations tend to grow 
faster than the sum fluctuation, leading towards a limit where |e| ^ |e|, or 5p ~ —5p. This implies tendency towards 
anticorrelation of matter and TM overdensities, or to segregation of the two sectors. This applies not only when the 
perturbations are weak: the fact that distinct matter and TM bodies repel each other for low accelerations implies 
that the process of segregation continues during the stages where the overdensities are high, (v) There is an added 
enhancement, due to the fact that once the fluctuations in the two sectors are anticorrelated, with maxima of each 
residing in minima of the other, each overdensity grows not only by the pull of its own gravity, but also by shepherding 
by the repulsive action of the flanking overdensities in the other sector, (vi) There can be strong coupling between 
the perturbations on different scales due to the nonlinearity inherent in MOND. 

There is much that remains to be studied, even within the narrow paradigm of the MTMS universe. An even richer 
prospect opens if one considers more general cases such as matter- TM asymmetric BIMOND versions, or asymmetric 
universes in which the two sectors are not equally populated, or do not have exactly the same properties. 

For example, we do not know why a universe that presumably started with a vanishing baryon-number density 
has ended up having more baryons than antibaryons, with the difference resulting in the present-day baryon density. 
This baryogenesis, as it is called, which requires a breakdown of symmetries (CP and baryon-number conservation), 
is thought to have occurred during some early phase of thermal disequilibrium. If the corresponding conditions in 
the two sectors where slightly different, we could have ended up with different baryon densities in the two sectors. It 



Even if the initial fluctuations arc different in the two sectors, we saw that their statistical properties are driven to equalization by the 
faster, MOND-like growth of the fluctuation difference. 
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would be interesting to study the effects of such departure from MTMS, to lowest order in the density difference, to 
see what effects it has on both the cosmological evolution, and on structure growth. 
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